Effective Field Theory for Goldstone Bosons in Nonrelativistic Superfluids 
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We consider nonrelativistic superfluids where the global I7(l)-symmetry is spontaneously broken. 
At sufficiently long wavelengths, the relevant degree of freedom is the massless Goldstone mode and 
we construct an effective low energy theory for the Goldstone boson. The damping rate of collective 
excitations at low energy is calculated. In the case of a weakly interacting Bose gas, we recover the 
results by Beliaev, and by Hohenberg and Martin. 
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I. INTRODUCTION 

Spontaneous symmetry breaking is central in our de- 
scription of many phenomena in condensed matter and 
high-energy physics (see e.g. Ref. [|| and references 
therein). An example of spontaneous symmetry break- 
ing is the breaking of the global [/(l)-symmetry in ^He 
and the appearance of superfluidity below a certain crit- 
ical temperature. Associated with the breaking of the 
symmetry, there is a nonzero value of an order parame- 
ter and a condensate of particles in the zero-momentum 
state. When a global continuous symmetry is broken, the 
Goldstone theorem |^ states that there is a gapless exci- 
tation for each generator that does not leave the ground 
state invariant. In the case of nonrelativistic Bose gases, 
one identifies the Goldstone mode with the phonons. 

At long wavelengths where the relevant degrees of 
freedom are the massless Goldstone bosons, one would 
naturally like to have an effective low energy theory that 
contains only these degrees of freedom. Such an effec- 
tive field theory may simplify calculations of the long- 
wavelength properties of the system and one can obtain 
model independent predictions with a minimum of as- 
sumptions. 

An example of a low energy theory for mass- 
less Goldstone modes is chiral perturbation theory in 
QCD 1^. In massless QCD, chiral symmetry is broken 
and one identifies the Goldstone bosons with the pions. 
However, QCD is a confining and strongly interacting 
theory at low energies. Thus the coefficients of the chiral 
Lagrangian cannot be calculated from QCD. Instead, the 
coefficients of the effective Lagrangian are fixed by exper- 
iments. In other cases, one can determine the coefficients 
of the low-energy theory as functions of the coupling con- 
stants in the underlying theory e.g. by a matching proce- 
dure. One calculates physical quantities at low energies 
perturbatively and demand they be the same in the full 
and in the effective theory. The coefficients of the effec- 



five theory then encode the physics at short distances. 
Nonrelativistic QED is an example of such an effec- 
tive field theory that is taylored to perform low-energy 
(bound state) calculations, where the coefficients encode 
the effects of relativistic momenta. 

In this paper, we construct an effective low energy 
theory for the Goldstone boson in nonrelativistic super- 
fluids with a broken global L/(l)-symnietry. The effective 
Lagrangian allows one to calculate physical quantities in 
the long-wavelength limit with ease without assuming 
weak coupling. 



II. EFFECTIVE LAGRANGIAN 

The effective theory for the Goldstone field (f) can be 
constructed using the methods of effective field theory [|j . 
Once the symmetries of the theory have been identified, 
one writes down the most general Lagrangian C^s that is 
consistent with these symmetries. Examples of symme- 
tries are time-reversal invariance and Galilean invariance, 
and these symmetries severely restrict the possible terms 
in £eff- For instance, Galilean invariance implies that the 
terms in the effective theory are powers of the combina- 
tion idr(j> — \ (V'A)^ {h — kB — m — 1 henceforth). The 
Euclidean Lagrangian can then be written as 
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where Ci,C2, ... are parameters that must be determined 
either by experiment or by matching. The first term is 
a total derivative and can be omitted if one is not inter- 
ested in topologically nontrivial field configuration such 
as vortices. The term SCcs contains all terms that are 
higher order in the field (j) and derivatives thereof. We 
have also omitted the unit operator / which is nec- 
essary to include if one is interested in the equation of 
state. The unit operator / can be interpreted as the con- 
tribution to the free energy from large momenta. We will 



briefly discuss this issue below. The terms in Eq. that 
we have shown explicitly, gives rise to a linear dispersion 
relation for the Goldstone bosons, and the effective La- 
grangian is therefore valid only for momenta where this 
is a good approximation. At larger momenta, the dis- 
persion relation is no longer linear and terms of the form 
(V^(/))^ must be included. The scale at which the effec- 
tive Lagrangian no longer can be applied is given by the 
scale where other degrees of freedom than the Goldstone 
boson become important. In the case of a weakly in- 
teracting Bose gas, this scale is given by the coherence 
length ^. 

The free propagator that corresponds to the La- 
grangian (0) is 



A(po,p) 



1 
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(2) 



where pa = 27rnT are the Matsubara frequencies and 
p = |p|. The dispersion relation is given by the pole of 
the Minkowski space propagator: 



e(p) = \I-P- 



(3) 



The nonrelativi stic h ydrodynamic speed of sound c is 
then given by ^ cijcx. In some cases, one does not know 
the underlying microscopic theory, and so one cannot de- 
termine the couplings ci, C2, ... by matching. Instead, one 
performs experiments at low energies. For instance, the 
ratio C2/C1 is fixed by measuring the speed of sound. 

In the case of a weakly interacting Bose gas, we know 
the underlying theory and we can calculate the coeffi- 
cients ci, C2, ... in terms of the parameters of the full the- 
ory. The Euclidean action for nonrelativistic bosons at 
low energy is 



d_ 



(4) 



where ipi'^) annihilates a boson at position x, /i is 
the chemical potential, and g — Ana, where a is the 
s-wave scattering length. If we substitute t) = 
y^TiQ + cr(x, T)e'^'^^^''^^ into Eq. (||), we obtain the action 



S[a,<j)] = J dr J dxUnog - fi]nocr + 



dcj) , da 


OT OT 



1 2 1 

-2"o(V(/>) +-(T 



— V" + Vo 

4no 



(5) 



where we have dropped total derivatives as well as an 
infinite series of higher order momentum-dependent in- 
teractions. By using the classical equation of motion for 
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FIG. 1: One-loop Feynman diagrams for the self-energy 

n(po,p). 



(t(x, r), we obtain from Eq. 
the phase 4>: 



the following action for 



S[(t)] = J^dT J dx 



.(6) 



Using Eq. ([|) and S = Jq dr J dx CcB, and comparing 
with Eq. (|) we find the coefRcients ci = no and C2 — l/g. 

As an application of the effective low-energy La- 
grangian (0) for the Goldstone bosons, we next calculate 
the damping rate of a long-wavelength excitation in the 
superfluid. The propagator has poles in the complex en- 
ergy plane, and the real part gives the energy e{p) of the 
excitation, while the imaginary part gives the damping 
rate. 

The effective Lagrangian gives rise to three-point 
and four-point interactions and the one-loop Feynman di- 
agrams contributing to the self-energy Tl{po, p) are shown 
in Fig. |l|. The first diagram vanishes identically since the 
vertices are momentum dependent. The second and third 
diagrams are ultraviolet divergent and must be regular- 
ized. In this paper, we use dimensional regularization 
to regulate both infrared and ultraviolet divergences. In 
dimensional regularization, one calculates the loop inte- 
grals in d = 3 — 2e dimensions for values of e where the 
integrals converge. One then analytically continues back 
to d = 3 dimensions. With dimensional regularization, 
an arbitrary renormalization scale M is introduced. An 
advantage of dimensional regularization is that it auto- 
matically sets power divergences to zero, while logarith- 
mic divergences show up as poles in e. The second di- 
agram is proportional to p^ and has pentic, cubic, and 
linear divergences in the ultraviolet. It is therefore set 
to zero with dimensional regularization. The third graph 
has a logarithmic ultraviolet divergence. The one-loop 
self-energy is 
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where 7 « 0.5772 is the Euler-Mascharoni constant. In 



2 



the zero-temperature hmit, the sum over Matsubara fre- 
quencies becomes an integral over the Euclidean energy. 
The next step in evaluating the self-energy is to introduce 
a Feynman parameter y. After integrating over energy, 
momentum, and finally over y, the self-energy can be 
expanded in powers of e. The result is 



n(po,p) = - 
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where g{po,p) is an analytic function that is not im- 
portant to us. After analytic continuation to real fre- 
quencies, po ^ iuj + 77, we find the the damping rate 
7 = Imn(-ie(p) + rj, p)/2e(p) 
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(9) 



Note that the result is independent of the coupling C2. By 
measuring the damping of long- wavelength excitation one 
can determine ci, while measuring the speed of sound c 
will determine ratio cxjci. Using the value c\ — hq for 
a weakly interacting Bose gas, we recover the result by 
Beliaev § (See also Refs. |, 

The quadractic part of action has recently been 
used a starting point for calculating phase fluctuations 
in trapped Bose gases at low temperature |10, O. The 



effective field theory approach clearly shows why such 
an approach is incorrect. The phase-alone action is only 
valid for momenta less than the inverse coherence length, 
while in Refs. 0, it was used for all momenta. The 
correct treatment of the phase fluctuations was given in 
Refs. @||. 

Consider next the pressure for a dilute Bose gas at 
zero temperature. If one calculates V using the effec- 
tive Lagrangian one gets zero. This disagrees with 
the one-loop result of Lee and Yang Q. The reason 
is simply that one must include the unit operator / in 
the effective Lagrangian (Q) and that contribution to the 
pressure exactly equals the standard result. 



III. FINITE TEMPERATURE EFFECTS 

We next consider finite-temperature effects using the 
effective low-energy Lagrangian (|^) 

The pressure at one-loop order is given by 
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(10) 



For a weakly interacting Bose gas, this result is in agree- 
ment with the leading temperature-dependent term in 



a low-temperature expansion carried out by Lee and 
Yang |l| (See also Ref. This simply reflects the 

fact that the massless Goldstone modes are the impor- 
tant thermal excitations at low temperature. 

Finally we consider damping at finite T. The imagi- 
nary part of the self-energy can be calculated in a similar 
manner as before, except that the integral over Euclidean 
energy is replaced by a summation over Matsubara fre- 
quencies. The result is 
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In the case of a weakly interacting Bose gas, this results 
was first obtained by Hohenberg and Martin 1 17 . Using 
the Lagrangian (|), Liu |] obtained a more general result 
to lowest order in p. This result is valid for p <^ T 
and p <C n^g and reduces to iii^ pT'^ / [AQno) in the limit 
T/nog 0. 



IV. DISCUSSION 

In this paper, we have constructed an effective field 
theory for Goldstone modes at low energies for a nonrela- 
tivistic superfluid with a broken [/(l)-symmetry. The re- 
sults for the damping of collective excitations is in agree- 
ment with those of Beliaev, and by Hohenberg and Mar- 
tin for the weakly interacting theory. However, the re- 
sults go beyond because nowhere have we assumed weak 
coupling. The only assumptions are Galilean invariance 
and low energy. Hence, the results (^ and ( pT| ) are valid 
for any nonrelativistic superfiuid with a broken U{1)- 
symmetry in the long- wavelength limit. 

There is another way of arrriving at a low-energy 
theory for the Goldstone modes, namely by calculating 
the quantum effective action ||l|. In the effective ac- 
tion approach, integrating out "short-distance physics" 
amounts to minimizing with respect to the amplitude of 
the field ^/'(x, r). Such an approach was used in Ref. Q 
to calculate a low-energy effective Lagrangian for rela- 
tivistic superfluids where the f/(l) baryon symmetry is 
spontaneously broken. In Ref. it was shown that 
the scattering amplitudes among the Goldstone bosons 
can be found once the equation of state is known. Ap- 
plying the formalism to the dilute Bose calculating the 
quantum effective action and using 7'(/i) — ^j? jlg^ one 
would immediately obtain Eq. (|l|) with the coefficients 
ci = /.t/g = TT-o and C2 = l/g. 
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